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Abstract. 

Within a first principles framework, we study the electronic structure of the 
recently synthesized polymeric coordination compound Cu(II)-2,5-bis(pyrazol-l-yl)- 
1,4-dihydroxybenzene (CuCCP), which has been suggested to be a good reaHzation of 
a Heisenberg spin- 1/2 chain with antiferromagnetic coupHng. By using a combination 
of classical with ab initio quantum mechanical methods, we design on the computer 
reliable modified structures of CuCCP aimed at studying effects of Cu-Cu coupHng 
strength variations on this spin-1/2 system. For this purpose, we performed two types 
of modifications on CuCCP. In one case, we replaced H in the linker by i) an electron 
donating group (NH2) and ii) an electron withdrawing group (CN), while the other 
modification consisted in adding H2O and NH3 molecules in the structure which change 
the local coordination of the Cu(II) ions. With the NMTO-downfolding method we 
provide a quantitative analysis of the modified electronic structure and the nature of 
the Cu-Cu interaction paths in these new structures and discuss its impHcations for 
the underlying microscopic model. 
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1. Introduction 

Low-dimensional quantum spin systems have been a subject of intensive research in 
the last decades, especially in solid state theory. The discovery in recent years of low- 
dimensional materials has allowed a direct comparison of experimental observations with 
theory predictions. Out of these studies, the need of available low-dimensional materials 
with a certain degree of tunability has emerged since tunability has the advantage of 
allowing a step by step understanding of the complexity of the materials. In this paper 
we will focus on the search and study of tunable low-dimensional materials. We will 
restrict ourselves to a particular class of systems, namely metalorganic compounds. 

Metalorganic compounds formed by transition metal centres bridged with organic 
ligands are being intensively discussed in the context of new molecule-based magnets 
and electronic materials [H [2]. They constitute a class of tunable materials partly 
due to their modular nature. The modular setup has the advantage of allowing the 
modification of relevant subunits chemically without changing the subsequent crystal 
engineering. Substitution of organic groups and ligands in these systems play the role 
of doping in the search for materials with desired magnetic interaction strengths and 
charge carrier concentrations. 

In the present work we pursue these ideas from a theoretical point of view. We 
consider a computationally feasible combination of classical with quantum mechanical 
ab initio tools [3J in order to design and analyze new metalorganic compounds. As an 
example, we introduce systematic changes on existing metalorganic materials in order to 
achieve desirable electronic or magnetic properties in the modified new structures. Such 
study i) allows for a gradual understanding of the properties of these low-dimensional 
systems ii) provides a guide to systematic synthesis in the laboratory. 

We focus our attention on the recently synthesized coordination polymer Cu(II)- 
2,5-bis(pyrazol-l-yl)-l,4-dihydroxybenzene [4] (CuCCP) which, from susceptibility 
measurements [5], has been identified as a model system for a spin- 1/2 Heisenberg chain 
with an antiferromagnetic exchange coupling constant of J = 21.5 K. The polymeric 
unit is shown in Fig. [U This compound is a good starting point to study effects of 
coupling strength variation by appropriately introducing modifications on the linkers 
as well as changes on the Cu coordination. Note that the coupling strength scale for 
this system is small compared f.i. to inorganic Cu oxides [6], where the couplings are 
more than one order of magnitude larger. Accordingly, in this work we will be expecting 
coupling strength variations in the range of one to a few tens of meV in energy. 

The Cu-Cu interaction in this compound depends on the electronic nature of 
the linker. Its properties can be tuned smoothly and predictably by changing the 
substitution pattern R (see Fig. [T]) of the central benzene ring (hydroquinone), or by 
introduction of additional ligands X at the Cu(II) ions. The substitution or introduction 
of additional ligands is expected to bring changes in the electronic properties of the 
compound. For example, it can change the magnitude of magnetic interactions between 
the Cu(II) centres in the spin chain via change in the charge density in the polymeric 
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Figure 1. Polymeric unit of Cu(II)-2,5-bis(pyrazol-l-yl)-l,4-dihydroxybenzene 
(CuCCP) (X=0, R=H). We will consider the substitutions R=CN and R=NH2 on 
the central benzene ring and the ligands X=H2 and X=NH3. 

chain. It may change the effective, inter-chain interactions, the one-dimensional nature 
of the original compound may thus be modified. It may even change the coordination 
and valence of the Cu(II) ions, which may induce changes in the transport properties 
along the one-dimensional chain by moving away from the Mott insulator at a half-filled 
Cu ?>dx'2-y-2 orbital. 

With the above ideas in mind, we first considered two possible H replacements 
in the central benzene ring: an NH2 group, which acts as electron donating group, 
and a CN group, which acts as electron withdrawing group (see Fig. [T]). Furthermore, 
extensive crystallization trials showed that crystallites of the Cu(II) polymer always 
contain lattice defects in high concentration. In many compounds, the Cu(II) ions are 
coordinated by six nitrogen or oxygen atoms instead of four ligands. The Cu(II) polymer 
is experimentally crystallized in a mixture of water and ammonia solvents, and it is likely 
that H2O or NH3 molecules are built into the crystal lattice. We have, therefore also 
considered the introduction of additional ligands like H2O or NH3 at the Cu(II) ions 
in our simulation study. In Fig. [T] and Table [T] we give an account of these modified 
structures. 

The paper is organized as follows: In section [2] we give an account of the 
computational methods considered in this work. In section [3] we present our results 
and discussion of the relaxed new crystal structures and the analysis of their electronic 
properties. Finally in section [4] we summarize our findings. 
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Short name 


R 


X 


CuCCP 


H 


- 


Cu(II)-NH2 


NH2 




Cu(II)-CN 


CN 




Cu(II)-H20 


H 


H2O 


Cu(II)-NH3 


H 


NH3 



Table 1. Naming convention for the substitutions and ligands on the CuCCP 
coordination polymer considered in this work. 

2. Method 

The methods used in this work, can be primarily categorized into two classes. Firstly, 
a class of methods has been used for the accurate structural determination of both the 
parent and the modified compounds. Once the structural aspects are decided, their 
electronic structures are calculated and analyzed with another class of methods. Note 
that the understanding of a complex system and design of new compounds need a 
combination of several different methods, each being focused to deal with one specific 
aspect. In the following, we give a brief description of all the methods that we have 
employed. 

2.1. Determination of Crystal Structure 

In the absence of diffraction data, a method much used to a priori predict crystal 
structures, is the force field technique [Zl El El [IDl ttH [E]- While such calculations 
are computationally fast, they rely on a classical ansatz and therefore miss all possible 
quantum mechanical effects, which are important for the description of the electronic 
structure. Quantum mechanical methods, on the other side, are computationally 
much more demanding, and they are, in this context, typically employed for two 
tasks: One is the local optimization after global optimization with force field methods. 
This has been reported for inorganic systems like NaCl or MgF2 [El [E] and for 
simple organic compounds like glycole C2H4(OH)2 and glycerol C3H5(OH)3 [T5] . 
Another is for secondary computations like the determination of molecular geometries, 
electrostatic charges or for the calculation of intramolecular and intermolecular potential 
curves [El [E] • 

In the present work we used an effective way of designing reliable crystal structures 
which shares the advantages of both methods, namely the fast calculations with classical 
force field methods and the subsequent accurate quantum mechanical description with ah 
initio methods. We first created the modified structures on the basis of crystallographic 
databases [E] and the crystal structures were optimized by force field methods. In 
the second step the structures were relaxed by ah initio quantum mechanical molecular 
dynamics |19) within the density functional theory formalism (DFT) until the forces on 
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the atoms were less than a given threshold to ensure structure stability. Our work differs 
from the mentioned previous works in the sense that using this approach, we succeeded 
in treating materials with large unit cells (of the order of 100 atoms) and complicated 
electronic structure (transition metal complexes) with sufficient accuracy. 

All force field optimizations were performed using the program package Cerius2 [20] . 
We modified the Dreiding 2.21 plj force field by introducing energy terms for the case of 
octahedrally coordinated metal ions. For the energy minimizations we used the modified 
Dreiding force field with Gasteiger [22] charges. All structural models were based on 
the experimentally determined crystal structure of CuCCP [4]. The crystallographic 
symmetry of the structure models was maintained in all relaxations. The position of 
the Cu ion was kept fixed during all force field and quantum mechanical optimizations. 
The second step of quantum mechanical relaxations were performed by Car Parrinello 
(CP) ah initio molecular dynamics (AIMD) calculations [T9] based on the Projector 
Augmented Wave (PAW) method [231124] . 

2.2. Electronic Structure Calculations 

We computed the electronic structure of the relaxed structures with the Full Potential 
Linearized Augmented Plane Wave basis (FPLAPW) as implemented in the Wien2k 
code [25]. Calculations were performed within the Generalized Gradient Approximation 
(GGA) [26]. The choice of muffin-tin radii tmt, k mesh and plane- wave cutoffs /Cmax 
were carefully tested. We considered a k mesh of (8 x 6 x 5) in the irreducible Brillouin 
zone and a -Rfcmax = 3.8, which is reasonable for systems that contain hydrogen atoms. 

It is well known that LDA or GGA fails to describe the correct insulating ground 
state for strongly correlated electron system, as is the case here. Introduction of missing 
correlation effects in a static mean-field like treatment as is done in the so called LDA+U 
approach[27l [28], should give rise to the correct insulating state, as is supported by 
our calculations (not shown here). In the present paper we are mainly interested 
in estimating the effective one-electron hopping interactions which are well described 
within LDA or GGA. In fact, the use of DFT calculations to understand the chemistry 
of correlated materials is a well established method [29] . 

Finally, in order to analyze the computed electronic structure and to extract 
an effective microscopic Hamiltonian, we derived quantitatively the Cu-Cu hopping 
integrals within the NMTO- downfolding technique [30l [31] as implemented in the 
Stuttgart code [32l|33]. The downfolding technique provides a useful way to derive few- 
orbital Hamiltonians starting from a complicated full LDA Hamiltonian by integrating 
out degrees of freedom which are not of interest. This procedure naturally takes 
into account the renormalization effect due to the integrated-out orbitals by defining 
energy-selected, effective orbitals which serve as Wannier-like orbitals for the few- 
orbital Hamiltonian in the downfolded representation. The method provides a first- 
principles way for deriving a few-band, tight-binding Hamiltonian of the form Htb = 
J2ijtij{CiCj + h.c.) for a complex system, where the tij define the effective hopping 
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between the downfolded orbitals and c| (q) are electron creation (annihilation) operators 
on site i. This method has proved to be extremely successful for systems such as high- 
Tc cuprates [34], double perovskites [35] or low-dimensional inorganic quantum spin 
systems [Ml EH EHl EH IQ] . 

Such estimates of the effective hopping integrals are useful in defining the underlying 
low-energy magnetic model. More precisely, the one-electron effective Cu-Cu hopping 
integral, t can be related to the Cu-Cu magnetic exchange coupling interaction J via 
a second-order perturbative treatment within the framework a many-body Hubbard- 
like model. Assuming that these couplings are antiferromagnetic and neglecting 
ferromagnetic contributions, J can be estimated as Jafm ~ 4t^/f4fT where Ues is the 
effective onsite Coulomb repulsion on the Cu site. 

3. Results 

3.1. Crystal Structure 

Dinnebier et al. [4] reported the synthesis, and crystal structure determination 
of the CuCCP system obtained by layering a solution of 2,5-bis(pyrazol-l-yl)-l,4- 
dihydroxybenzene in CH2CI2 with a solution of CuBr2 in concentrated aqueous 
ammonia. The system crystallizes in the triclinic space group PI with 27 atoms per 
unit cell. This compound tends to form independent polymeric chains consisting of 
deprotonated 2,5-bis(pyrazol-l-yl)-l,4-dihydroxybenzene molecules bridged by Cu(II) 
ions with a 3c?^ configuration, which corresponds to a local spin 1/2. As shown in Fig. 
[2] the chain axes are oriented along the c axis of the crystal and the copper ions are 
located at (1/2, 1/2, 1/2), which is a centre of symmetry of the space group PI. The 
Cu-Cu distance along the approximate a axis is about 5.2 A while along the other two 
axes it is close to 8 A. These Cu(II) ions are coordinated in an almost square planar 
fashion by two (pyrazolyl) nitrogen atoms and two oxygen atoms of the deprotonated 
dihydroxybenzene groups. 

The crystal structure was determined from X-ray powder diffraction data; 
consequently the overall structure and the arrangement of the chains are reliable, but 
the individual atomic positions had an accuracy of only about 0.3 A. A DFT analysis 
of the forces p5] between the atoms shows that the experimentally determined CuCCP 
structure is still very unstable with forces of the order of 100 mRyd/oe or more for some 
atoms. We have therefore relaxed the atomic positions keeping the Cu position fixed 
with the AIMD method described in section 12. 1[ 

The Cu(II)-NH2 polymer and the Cu(II)-CN polymer were generated from the 
CuCCP polymer, by substituting the two hydrogen atoms of the benzene rings by amino 
(NH2) or cyano (CN) groups, respectively (see Fig. [T] and Table [I]). The Cu(II)-H20 
(Cu(II)-NH3) polymer was constructed from the CuCCP polymer by adding two water 
molecules (ammonia molecules) as additional ligands to the Cu(II) ion (see Figs. [Hand 
and Tabled]). In the original crystal structure the chains are quite densely stacked. 
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Figure 2. Crystal structure of the CuCCP polymer in two different orientations. The 
unit cell is shown in the figures (vectors a, 5, and c) . Note the arrangements of the Cu 
chains along the c direction. The various Cu-Cu interaction paths ti have been also 
drawn where the index i = 1, 2, 3, 7, 8 denotes the i-th nearest neighbour. 
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Table 2. Lattice parameters for the experimental crystal structure CuCCP and models 
Cu(II)-H20 and Cu(II)-NH3, optimized with force field methods. 



Structure 


a(A) 


6(A) 


c(A) 


an 




in 




CuCCP exp 


5.172 


7.959 


8.230 


118.221 


91.520 


100.148 


291.47 


Cu(II)-H20 


5.234 


11.249 


8.072 


117.611 


68.822 


127.155 


330.43 


Cu(II)-NH3 


5.459 


11.597 


8.349 


118.423 


68.840 


130.883 


350.49 



The introduction of the H2O (or NH3) molecules would either lead to unrealistically 
short contacts to the neighbouring chains, or to a considerable increase of the distances 
between the chains, resulting in an unrealistically loosely packed structure. Therefore 
the crystal structures of the Cu(II)-H20 and Cu(II)-NH3 polymers were fully optimized, 
including an optimization of the lattice parameters. Moreover, in order to achieve a 
better packing of the Cu(II) polymer chains with a favourable lattice energy, the Cu(II) 
chains shifted in the optimization process both sidewards as well as along the chain 
direction with respect to each other. The resulting unit cell parameters are shown in 
Table El 

All the modified structures were relaxed in the second step with the AIMD method 
until the forces on the atoms were sufficiently small to ensure stability of the quantum 
mechanical calculations. In the Appendix we present the relaxed structural data of all 
the modified structures. 

3.2. Electronic Structure and Effective Cu-Cu Interactions 

3.2.1. CuCCP In Fig. [3] we present the band structure for the relaxed CuCCP where 
the Cu d band character is shown by bars on the right side. The band characters are 
given in the local coordinate frame of Cu which is defined with the local z direction 
pointing from the Cu to out-of-plane N atom in the next layer and the y direction 
pointing from the Cu to in-plane O atom. Cu is in a "id^ configuration, with all d bands 
occupied except for the last band which is half-filled. GGA predicts a metallic behaviour 
for this system. As mentioned previously, inclusion of on-site electronic correlation 
within LDA+U opens a gap between a lower occupied Hubbard band and an upper 
unoccupied Hubbard band and the system is described as a Mott-Hubbard insulator. 
Since the 0-Cu-N angle in the CUO2N2 plane is not exactly 90°, the various Cu d 
degrees of freedom defined with respect to the local coordinate frame mentioned above 
show slight admixtures. In particular, the Cu dx2_y2 dominated band crossing the Fermi 
level has also small contributions from Cu dyz degrees of freedom which arise from the 
distorted geometry. 

>From the dispersion of the Cu d band at the Fermi level we confirm the one- 
dimensionality of the structure. The paths F-F and B-F which correspond to the inter- 
chain paths are almost dispersionless and the intrachain F-Z-B path shows a cosine-like 
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F r z Br 

Figure 3. Band structure for the relaxed Cu(II) polymer CuCCP in the 
GGA approximation along the path [41] F(0, 1, 0)-r(0, 0, 0)-Z(0, 0, 1)-B(0.99, -0.13, 0)- 
r(0, 0,0) in units of 7r/a, 7r/6, tt/c. The bars indicate the dominant band character in 
the local coordinate frame of Cu (see text for explanation) . 

behaviour. A quantitative analysis of the various hopping integrals obtained with the 
downfolding procedure, by keeping only the Cu dx2-y2 degrees of freedom active and 
integrating out all the rest within the NMTO framework, is given in the first column 
of Table [3l The various interaction paths are as shown in Fig. El The largest hopping 
integral is along the chain (see Fig. [2j) while all other hoppings are almost an order 
of magnitude smaller. 

3.2.2. Cu(II)-NH2 and Cu(II)-CN In order to have a quantitative account of the 
structural changes that the polymer system undergoes under the various substitutions, 
we define the angle between the vector perpendicular to the CUO2N2 plane and the 
vector perpendicular to the benzene ring as the tilting angle {}. The substitution of H 
by NH2 groups or CN groups in the benzene rings induces a tilting from d = 34.9° in 
CuCCP to^ = 37.3° in Cu(II)-NH2 and 1^ = 36.3° in Cu(II)-CN. 

In Figs. [1(a) and (b) we present the FPLAPW orbital resolved DOS for the Cu(II)- 
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Figure 4. Orbital resolved DOS for (a) Cu(II)-NH2 and (b) Cu(II)-CN. Panel (c) 
shows the comparison of the contribution of the NH2 and CN groups to the DOS at 
Ey in a blown up scale. 
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Table 3. Values for the Cu-Cu hopping integrals calculated with the NMTO 
downfolding method for the relaxed CuCCP, Cu(II)-NH2 and Cu(II)-CN structures. 
The values are given in meV. The subscripts i= 1,2,3,7,8,12 denote the ?th nearest 
neighbours. See Fig. [21 Only the hopping integrals having values larger or equal to 
one tenth of a meV has been shown. 
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NH2 and the Cu(II)-CN within the GGA approximation. Shown is the contribution to 
the total DOS of Cu, O, N, C and the groups NH2 and CN. While some changes in 
the detail shape of the DOS for Cu, O, N and C between Figs. [H (a) and [H (b) are 
observed, the most important effect is the different electronic nature of the NH2 and 
CN groups. Fig. [4] (c) shows the contribution of these groups at the Fermi surface. 
The CN group bands are deep down into the valence band while the NH2 group has 
appreciable contribution near the Fermi level, which indicates its involvement in the 
effective interaction paths between copper atoms. We will see this more clearly in the 
plot of the NMTO-Wannier orbitals to be discussed later in this section. 

In Fig. Owe show a comparison of the band structure for the relaxed CuCCP, Cu(II)- 
NH2 and the Cu(II)-CN polymers in the energy range [—0.25 eV, 0.25 eV] where only 
the Cu dx2-y2 dominated band is involved. Though the basic nature of the dispersion 
remains the same upon substitution, the details however do change. NH2 seems to be the 
most effective substitution to increase the intrachain Cu-Cu interaction (the bandwidth 
widens along the F-Z-B path for the Cu(II)-NH2 system). The CN substitution, on 
the other hand, reduces this interaction (note the bandwidth narrowing along the F- 
Z-B path for the Cu(II)-CN system). The substitution process also enhances certain 
interchain couplings. The almost dispersionless behaviour along F-F and B-F becomes 
more dispersive. Description of such fine and subtle changes, need some quantitative 
measures which can be best described by the changes in effective Cu-Cu hoppings. This 
is shown in Table [3] where the hopping integrals obtained by the NMTO downfolding 
method are shown. Note that the ti hopping for the Cu(II)-CN system along the 
crystallographic a direction is enhanced by a factor of 4.5. Similarly, tj and ts hopping 
terms for the Cu(II)-CN system between neighbouring Cu chains in the b direction (see 
Fig. [2]) are almost 2-3 times larger compared to that of the CuCCP system. The long- 
ranged ti2 hopping parameter between neighbouring chains along the a axis also attains 
appreciable enhancement compared to a vanishing small value for the CuCCP system. 
Similarly, ti and tg hoppings are enhanced for the NH2 substitution by factors of about 
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~ 2-3. Among all the hoppings, only t2 shows the exception of being systematically 
decreased upon substitution. The predominant hopping, ts is enhanced in the Cu(II)- 
NH2 system and reduced in the Cu(II)-CN system as already predicted from bandwidth 
arguments. 

A very helpful tool to understand the origin of these changes is the plot of effective 
Cu Wannier orbitals. In Fig. [6] we show the Wannier orbitals for the three Cu(II) 
systems presented so far. The plotted Wannier orbitals are obtained with the NMTO 
downfolding technique as introduced in the section 12. 2[ 

The effective Cu Wannier orbital has the expected Cu dx2-y2 symmetry at the 
central Cu site while the tails sitting at other sites are shaped according to the 
symmetries of the various integrated out orbitals like the rest of Cu d orbitals, O p, 
N p or C p. Note that the Cu d, O p and N p antibonding orbitals in the basic CUN2O2 
square plaquette remain similar in all three cases but the effective orbital distribution 
in the benzene ring is markedly different. The changes are most prominent for the NH2 
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Figure 6. Cu Wannier functions for (a) the relaxed CuCCP polymer, (b) the Cu(II)- 
NH2 polymer, and (c) the Cu(II)-CN polymer, (d) indicates the atom positions 
common to (a)-(c). The N-C-C-C-H chain of atoms appearing above the Wannier 
function belongs to the next layer. 



Cu-based metalorganic systems: an ah initio study of the electronic structure 14 

substituted case with the tails attaining appreciable weight at the sites in the benzene 
ring. We also notice the occurrence of weight at the NH2 assembly which is in accordance 
with the orbital resolved DOS study (see Fig. [4|). This leads to an enhancement in both 
intra- and some interchain Cu-Cu interactions, caused by the larger overlap of the 
effective orbitals. The enhancement happens via two different routes: one is due to the 
different tilting of the benzene ring compared to the original compound and the other 
one is the opening of additional interaction paths via the NH2 group which enhances the 
intrachain as well as interchain interactions ti and as can be seen in the quantitative 
estimates of the hopping interactions in Table [3l In the case of the CN substitution, 
the opening of an additional intrachain pathway is absent, which is reflected in the 
reduced intrachain (^3) hopping interaction. However the mechanism via the tilting of 
the benzene ring is still operative which is seen in the enhancement of several interchain 
coupling, specially ti. 

Relating the magnetic coupling interaction J with the effective hopping interaction 
t via a relationship Jafm ~ 4t^/[/efr, as discussed in Section 12.21 and choosing Ues 
to be 5 eV [42] we obtain the nearest neighbour coupling for CuCCP system to be 
<^AFM ~ 58 K which is somewhat larger than the experimental estimate [5] obtained 
by fitting susceptibility data to an effective nearest neighbour Heisenberg model, but 
remains of the same order of magnitude. The Jafm values estimated for Cu(II)-NH2 
and Cu(II)-CN systems are Jafm ~ 72 K and ~ 43 K respectively. 

3.2.3. Cu(II)-H20 and Cu(n)-NH^ In our second set of modifications we introduce 
two kinds of ligands, H2O and NH3 in the CuCCP system in the way presented in 
Section 13. 1[ Our goal is to study the effect of H2O and NH3 satellites on the CuCCP 
structure as well as to search for possible routes to change the Cu coordination from 
planar to octahedral. As explained in Section [3Tl in order to obtain realistic structures, 
the optimization with the force field method was done without keeping the original 
cell fixed, since that would force very short intermolecular distances between the H2O 
(NH3) moieties and the neighbouring chains. The force field optimized structures (see 
Table [2]) and subsequently relaxed with AIMD are characterized by Cu-0 (O of the 
H2O molecule) distances of dc^o = 2.17 A while the Cu-0 and Cu-N in-plane distances 
are dcno = 1-99 A and dc^N = 2.01 A, respectively. This corresponds to a distorted 
octahedron elongated along the CU-H2O direction. For the case of the NH3 ligands 
the Cu-N (N of the NH3 molecule) distances are dcuN = 2.14 A, while the Cu-0 and 
Cu-N in-plane distances are dcuo = 2.02 A and (icuN = 2.03 A, also giving rise to 
an elongated octahedron along the CU-NH3 direction. The ligands close to the Cu(II) 
centre also induce a tilting of the benzene rings with respect to the CUO2N2 plane. 
From the initial angle of ^ = 34.9° in CuCCP the tilting due to the H2O ligand is quite 
significant, leading to -d = 42.9° in Cu(II)-H20. The NH3 molecule, by contrast, leads 
to a lowering of this angle to 'd = 31.8° in Cu(II)-NH3. 

In order to quantify the effect of the H2O and NH3 ligands on the electronic 
properties of CuCCP, we show in Table [4] the values of the Cu-Cu hopping integrals 
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Table 4. Values for the Cu-Cu hopping integrals calculated with the NMTO 
downfolding method for the relaxed CuCCP, Cu(II)-H20 and Cu(II)-NH3 structures. 
The values are given in meV. The subscripts i — 1,2, 3, 7, 8, 12 denote the ith nearest 
neighbours. See Fig. [21 
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calculated with the NMTO downfolding method where the hopping parameters for the 
original CuCCP have been included for comparison. Note that the intrachain Cu-Cu 
coupling is reduced by a factor of 1.5-3.5 with the inclusion of both ligands. The 
reduction is especially significant with NH3. The only Cu-Cu interchain path that 
is enhanced is ti which is between Cu in nearest neighbour chains and has its origin 
in the hydrogen bonds between the H of the H2O (NH3) molecule and the O of the 
hydroquinone fragments in the chains. Therefore, apart from these hydrogen bonds, the 
introduction of ligands isolates the Cu ions considerably. 

In the Wannier orbital plot in Figs. [7] and [8] we can see that the distorted octahedral 
environment of the Cu in the Cu(II)-H20 and Cu(II)-NH3 structures induces very little 
mixing of the Cu orbital to the predominant d^^-y^- Also note the little contribution 
of weight in the hydroquinone ring, in contrast to the previous discussed systems (see 
Fig. [6]) which is a manifestation of the isolated nature of Cu in these structures. The 
inclusion of the H2O and NH3 satellites, however, changes the Cu coordination only 
marginally from four in the direction of six, as opposed to our original motivation for 
the addition of H2O and NH3 ligands. Cu remains in the oxidation state of almost 2+ 
as observed in our calculations. While the GGA calculations give a metallic behaviour 
with a half-filled predominantly dx2-y2 Cu band, inclusion of correlation effects with 
LDA+U drive the system to an insulating state. Therefore the system will remain an 
insulator. 

4. Summary 

In search of low-dimensional quantum spin systems with tunable properties, we 
have proposed and analyzed within an ab initio framework and using a combination 
of different computational methods various chemical modifications to the Cu-based 
polymeric coordination compound CuCCP. Our goal has been to tune in a controlled way 
the magnetic interactions between Cu centres and to test the efficiency and feasibility 
of the combination of methods proposed here. We pursued two ways of modifying 
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Figure 7. Cu Wannier functions for the Cu(II)-H20 system. 

the original CuCCP structure; by changing the substitution pattern in the linker 
(hydroquinone) and by adding ligands to the system. Following the first scheme we 
considered two possible H substitutions in the hydroquinone; an electron donating 
group (NH2) and an electron withdrawing group (CN). For the second scheme, we 
considered structures with H2O and NH3 ligands. Out of our study we conclude that 
the NH2 substitution in the hydroquinone is the most eff'ective in order to enhance the 
intrachain Cu-Cu interaction in CuCCP while the CN substitution induces and enhances 
the interchain interactions in the system which were either absent or very weak in the 
original CuCCP compound. In contrast, the inclusion of H2O or NH3 ligands has the 
effect of isolating the Cu ions. 

The effects observed in this study are small, mainly due to the fact that the coupling 
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Figure 8. Cu Wannier functions for the Cu(II)-NH3 system. 

constants in these metalorganic materials are weak. On the other hand, these systems 
are, due to the smallness of the coupling constants, of special interest since application of 
moderate magnetic fields or pressure can drive the system to a phase transition. These 
systems can be described as Mott-Hubbard insulators and possibly under application of 
pressure a metal-insulator transition could be induced. This will form the basis of our 
future work. 

Finally, we believe that the combination of methods presented in the manuscript is 
efficient for studying the properties of complex systems per se, and that the computer- 
designing procedure that we have employed in the present study provides a plausible 
route for manipulating properties related to low-dimensional quantum spin systems in 
general. 
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Table Al. Fractional atomic coordinates of nonequivalent atoms in the CuCCP 
relaxed structure. For the lattice parameters, see Table [21 



Atom 


X 


y 


z 


Cu 





5000 


0.5000 


0.5000 


02 





4581 


0.3495 


0.6187 


C3 





4745 


0.4243 


0.7998 


C4 





6696 


0.5869 


0.9293 


C5 





0331 


0.7729 


0.6797 


C6 





2206 


0.8455 


0.8363 


C7 





1024 


0.7795 


0.9497 


C8 





6971 


0.6615 


0.1203 


N9 





8576 


0.6739 


0.8610 


NIO 





8109 


0.6692 


0.6956 


Hll 





1732 


0.7923 


0.0815 


H12 





4193 


0.9317 


0.8607 


H13 





0456 


0.7875 


0.5549 


H14 





8471 


0.7919 


0.2121 



Appendix 

We present here the AIMD relaxed structural data of the various Cu(II) polymers. 
Other than CuCCP, the rest are computer designed. 

In Table lAll we show the fractional atomic positions of the CuCCP polymer 
obtained after the relaxation of the system. The resulting distances between the atoms 
after the optimization are closer to the standard values found in the literature [43] but 
differs from the values for the distances in the experimental compound [4]. 

In Table IA2I we present the relative atomic positions obtained after relaxation 
for Cu(II)-CN. The cell parameters were fixed during relaxation, and the only change 
observed was the tilting of the hydroquinone ring as a consequence of the movement 
upward of the pyrazolyl rings. 

In Table [A3] we show the same, but for Cu(II)-NH3. 

The structural data for the Cu(II)-NH2 and Cu(II)-H20 systems were already 
presented in Ref. [3j. 
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Table A2. Fractional atomic positions of nonequivalent atoms in Cu(II)-CN obtained 
with the PAW method. For the lattice parameters, see Table [2] (same as CuCCP). 



Atom 


X 


y 


z 


Cu 





5000 


0.5000 


0.5000 


02 





4559 


0.3460 


0.6195 


C3 





4823 


0.4252 


0.7998 


C4 





6795 


0.5881 


0.9253 


C5 





0384 


0.7593 


0.6620 


C6 





2349 


0.8314 


0.8139 


C7 





1239 


0.7694 


0.9310 


C8 





6944 


0.6601 


0.1182 


C9 





1030 


0.1811 


0.7598 


NIO 





8687 


0.6739 


0.8553 


Nil 





8152 


0.6676 


0.6880 


N12 





9162 


0.0575 


0.6788 


H13 





2104 


0.7793 


0.0570 


H14 





4353 


0.9135 


0.8323 


H15 





0493 


0.7681 


0.5346 
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